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1 INTRODUCTION 

Low-mass X-ray binaries (LMXBs) consist of a neutron 
star or a black hole in a pair with an optical star of mass M < 
IMq. The companion star in a close binary fills its Roche 
lobe and outflows onto a compact object with the formation 
of an accretion disk. 

The possibility of irradiation of the outer accretion disk 
in X-ray binaries was first pointed out by Shakura and Sun- 
yaev (1973). Since the outer accretion disk has a curved 
surface (H oc R", n = 9/8 for a standard disk), the fraction 
of the X-ray flux from the central source that is absorbed 
and "thermalized" near the disk surface depends on radius 
as 2* cx 1 /R-. At the same time, the intrinsic energy release 
(through viscosity) in the disk decreases with radius faster: 
oc For a binary with a neutron star, the disk sur- 

face heating through irradiation exceeds the intrinsic energy 
release already at radii ^ > 10^ cm. 

The vertical structure of the outer accretion disk in 
LMXBs with direct irradiation from the central source was 
investigated by many authors (see, e.g., Tuchman et al. 
1990; Vrtilek et al. 1990; Dubus et al. 1999). The X-ray 
flux is commonly assumed to be thermalized near the accre- 
tion disk photosphere (see, e.g., Dubus et al. 1999). In this 
case, external irradiation has no significant influence on the 
vertical structure of the disk as a whole (Lyutyi and Sunyaev 
1976). Indeed, surface heating affects the temperature in the 
central plane of the accretion disk if the following condition 
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is met: 



where tq is the total optical depth of the disk in the vertical 
direction for intrinsic radiation. The standard disk (Shakura 
and Sunyaev 1973) has a fairly large optical depth in the 
vertical direction (tq > 500). The surface heating of such a 
disk cannot have a significant effect on the conditions in the 
central plane up to very large radii. 

The model of a vertically isothermal accretion disk (Vr- 
tilek et al. 1990) is commonly used to describe LMXB ob- 
servations. For example, the disk thickness at the outer edge 
was estimated in this model from the observed amplitude of 
optical variations on the light curves of LMXBs (de Jong 
et al. 1996), H/R w 0.2, which is a factor of 2 ^ 3 larger 
than the thickness of an unirradiated standard disk (Shakura 
and Sunyaev 1973). Note that using the simple model of an 
isothermal disk (Vrtilek et al. 1990) to describe observations 
is inconsistent with the assumption about X-ray heating of 
only the surface. As has been noted above, such heating has 
no significant effect on the vertical disk structure, including 
the temperature profile. 

The necessity of an increase in the accretion disk thick- 
ness at the outer edge by a factor of 1.5-^2 compared to 
the standard disk thickness also follows from the model- 
ing of the observed X-ray light curves (Esin et al. 2000; 
Suleimanov et al. 2008). Note that this excess thickness can 
be explained by the presence of a hot atmosphere (see be- 
low) or a population of relatively cold clouds embedded in 
the atmosphere (Suleimanov et al. 2003) above the disk. 
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In an unirradiated accretion disk at a radius R ^ 5 x 
10"^ cm, the temperature on the photosphere drops below 
10000 K. As a result, a zone with incomplete hydrogen ion- 
ization appears in the disk and thermal instability must dis- 
rupt the stationary regime of accretion. On the other hand, 
as follows from observations, the semimajor axis of the or- 
bit in persistent long-period LMXBs with known parame- 
ters (see, e.g.. Table 1 in Gilfanov and Arefyev 2005) can 
be a w (2 ^ 5) X lO^^ cm (the binaries GX 13-h1 and Cyg X- 
2). The outer radius of the stationary accretion disk in such 
binaries. Rout ~ (2 8) x 10^' cm , is considerably larger 
than the radius at which a zone with incomplete hydrogen 
ionization appears in an unirradiated disk. 

The disk heating through X-ray irradiation, in principle, 
can increase the disk temperature, thereby moving the re- 
gion with incomplete hydrogen ionization farther along the 
disk radius. However, as was shown previously (see, e.g., 
Dubus et al. 1999), the central X-ray source cannot directly 
irradiate the outer disk region where hydrogen transformed 
into a neutral state, because the disk thickness in this zone 
decreases sharply (due to an increase in themolecular weight 
of the material during hydrogen recombination). The prob- 
lem with the screening of the outer accretion disk can be re- 
solved by assuming the presence of a scattering atmosphere 
with a temperature Tam « (2 3)10*' K above the disk. 

X-ray irradiation changes greatly the structure of the near- 
surface layers in an accretion disk. A thick hot (T ~ (10^ 
10^) K) plasma layer, which we call an atmosphere, can be 
formed above its surface. The calculations of a hot extended 
atmosphere appUcable to the accretion disks in LMXBs 
were performed in a number of works (Raymond 1993; Ko 
and Kallman 1994; Suleimanov et al. 1999; Jimenez-Garate 
et al. 2002). 

As was shown by Jimenez-Garate et al. (2002), there ex- 
ists a feedback mechanism between the disk irradiation and 
the thickness of its atmospheric layer, which leads to a sig- 
nificant increase of the latter. The atmospheric layer is op- 
tically thin in the vertical direction but optically thick in ra- 
dial coordinate. Thus, for X-ray photons from the central 
source, the thickness of the outer accretion disk increases by 
the height of the atmosphere. For example, in the above 
paper, the half-thickness of a disk with an atmosphere at 
R = lO'^ cm was found to be Zatm/R ~ 0.11 H-0.20 for an 
accretion rate in the range M = (0.1 ^ l)A^Edd> while the 
disk half-thickness at the photospheric level was Zph/R ^ 
0.062 ~ 0.083 in the same range of M. The presence of an 
extended atmosphere above the disk (as in Jimenez-Garate 
et al. 2002) naturally explains the enhanced disk thickness 
at the outer edge H/R ss 0.2 found by de Jong et al. (1996). 

Here, we constructed a self-similar model of an irradiated 
accretion disk by taking into account the scattering of X-ray 
photons in the disk. We also took into account the possibil- 
ity that an extended atmosphere was present above the disk 
(we took the atmospheric parameters from Jimenez-Garate 
et al. (2002), see §2.1). We expect that the Xray photons 



after their scattering in the disk and the atmosphere can pen- 
etrate and be thermahzed deeper in the disk than the direct 
photons from the central source. As will be shown below, 
the scattering effect leads to deep heating of all disk layers 
in z coordinate at a radius 7? ~ 10' ^ cm. In this case, the disk 
still remains optically thick for intrinsic radiation. 

The paper is organized as follows. In Section 2, we con- 
sider our model for calculating the vertical structure of an 
irradiated disk that includes the calculation of the transfer of 
X-ray photons in the disk and the atmosphere by the Sobolev 
method (Sect. 2.2) and the determination of the vertical disk 
structure under the photosphere (Sect. 2.3). The results of 
our calculation of the vertical structure of an irradiated ac- 
cretion disk and their discussion are presented in Section 3. 
In the final section, we summjirize our conclusions. 

2 CALCULATING THE VERTICAL STRUCTURE OF 
THE OUTER ACCRETION DISK WITH 
IRRADIATION FROM A CENTRAL X-RAY 

SOURCE 

2.1 Irradiation Parameters 

Consider a geometrically thin, optically thick accretion 
disk around a neutron star of mass M = I AMq in a LMXB. 
We assume the accretion rate in the disk to be constant: 
M = const. The chemical composition of the disk cor- 
responds to the cosmic elemental abundances: X=0.73, 
Y=0.25, Z=0.017 (Allen 1973). 

The X-ray radiation comes from a region near the com- 
pact object (a point source for the outer disk) and has an en- 
ergy spectrum S(E) corresponding to bremsstrahlung with a 
temperature Tsp = 8 keV. In the energy range 1 20 keV, the 
spectrum is 



-0.4 



expl-- 



(2) 



This model spectrum corresponds to the observed X-ray 
spectra of persistent LMXBs that can be described by 
a bremsstrahlung spectrum with a temperature T^p = 5 
10 keV (Liu et al. 2001). The disk-neutron star boundary 
layer is probably the main source of the radiation illuminat- 
ing the outer accretion disk regions. 
The X-ray luminosity of the central source is 



(3) 



where M is the mass accretion rate, the accretion efficiency 
onto the neutron star is assumed to be ~ 0.1. 

An X-ray flux (the central source is assumed to be a point- 
hke and isotropically emitting one). 



(4) 



is incident on the outer disk at radius R. Note that Eq. (4) 
is approximate. Thus, assuming the Xray radiation to come 
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Fig. 1. Irradiation geometry of a disk with an atmosphere and the coordi- 
nate system in the disk. 



from a boundary layer with a thickness much smaller than 
the neutron star radius without relativistic effects, 2it^ must 
appear in the denominator of Eq. (4) instead of the coeffi- 
cient Ait. 

The angle between the direction of incidence of the X-ray 
photons and the inward normal to the disk/atmospheric layer 
surface (see Fig. 1) is 

e,= ^--arctg{^-^+arctg{^y (5) 

where the function H{R) represents the dependence of the 
disk half-thickness (or the atmosphere height) on radius for 
X-ray photons. 

If the shape of the geometrically thin (H/R <^ 1) 
disk/atmosphere can be described by a power-law depen- 
dence on radius H(R}cxR", then the irradiation parameter 
Co = cos 6o is expressed as 



dH H H 



(6) 



We considered two models for the irradiation of the outer 
accretion disk: with and without an atmospheric layer above 
the disk. 

(1 ) The irradiation of a disk without an atmosphere. 

In this case, the disk thickness for X-ray photons corre- 
sponds to the photospheric level H = Zph for intrinsic radia- 
tion. The irradiation parameter Co depends on disk thickness 



and index nph (see Eq. (6)). We found zph simultaneously 
with the solution of the equations for the vertical disk struc- 
ture (see below). We fixed the power-law slope of the sur- 
face at riph = 9/8 (which corresponds to the outer parts of an 
unirradiated standard disk). 

(2) The irradiation of a disk with an atmosphere. 

There is an atmospheric layer that is optically thin in the 
vertical direction (r <C 1), but optically thick in the radial di- 
rection above the disk photosphere (for the possibility of the 
appearance of such a layer above the outer accretion disk, 
see Jimenez-Garate et al. 2002). 

A detailed calculation of the structure of an extended 
layer above an accretion disk is a nontrivial problem that 
is beyond the scope of this paper. To specify the parame- 
ters of the atmosphere above the disk, we used the results 
of atmospheric-layer calculations from Jimenez-Garate et 
al. (2002) obtained with similar irradiation parameters and 
chemical composition of the accretion disk. In the above pa- 
per, the atmosphere was assumed to be in hydrostatic equi- 
Ubrium. The stationary model atmosphere is applicable for 
most LMXBs with X-ray emission lines from a photoion- 
ized plasma observed in their spectra (X-ray pulsars and X- 
ray binaries with orbital inclinations i > 70°), except for bi- 
naries with a strong outflow of matter from the accretion 
disk (for more detail, see the discussion in §8.1 of Jimenez- 
Garate et al. (2002)). 

The disk thickness for X-ray photons corresponds to the 
photospheric level H = Zam- We assumed Zam and the slope 
of the surface of the atmospheric layer natm at a given radius 
to be (see Jimenez-Garate et al. 2001) 



nlR = 



1.2 X 10-2/?° 1^ 
l.Ox lO-^/jO^' 



,L/Lem ■ 
,L/LEdd = 



= 0.1 
1.0. 



(7) 



The total surface density in the atmospheric layer above the 
disk (see Fig. 14 in Jimenez-Garate et al. 2002) is 



0.066 
0.32 



■,L/Lem '■ 
,L/LEdd '■ 



0.1 

: 1.0 



(8) 



According to the calculations by Jimenez-Garate et al. 
(2002), the bulk of the surface density in the atmospheric 
layer is accumulated in a zone with a temperature T « 
2 3 X 10^ K, where the balance between heating and cool- 
ing is maintained by the photoionization/recombination of 
hydrogenand helium-Uke ions. A more rarefied region with 
a temperature T > 2 x 10' K, where the Compton heating 
and coohng processes dominate, Ues higher. Its contribution 
to Eatm is, on average, ~ 15% at radii R = 10'° 10" cm. 
As follows from the cited calculation, the total surface den- 
sity in the atmospheric layer T,atm is virtually independent of 
the radius in this range of R. 

In our simplified model of the atmospheric layer, we as- 
sumed its temperature to be Tam ~ 2 x 10^ K. The pressure 
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near the base of an isothermal atmosphere is 



where 3? is the universal gas constant, = \/ GMjB? is the 
angular Keplerian velocity at radius R, Ham is the molecular 
weight of the material, Zs is the height near the base of the 
atmospheric layer, zo = \J^^^^ and the function $(x) is 



2 pe-'^jr 

Jx 



Since the height near the base of the atmospheric layer w 
Zph, we obtain ^ ^ 1 and <&(Zj/zo) ~ 1- Thus, the formula 
for the pressure near the base of the atmosphere takes the 
form ^ 

'^^^ ^j, . (9) 

The molecular weight of the material in the atmospheric 
layer (a completely ionized medium) is /ia„„ = 0.6. 

2.2 Calculating the X-ray Radiation Transfer in the Disk 
and the Atmospheric Layer in the Sobolev 
Approximation 

Being absorbed and thermalized in subphotospheric lay- 
ers, the X-ray radiation can serve as an additional disk heat- 
ing source. 

To determine the mean intensity and flux of Xray photons 
in the accretion disk, we solved the onedimensional problem 
of X-ray photon transfer in a layer of finite thickness in the 
Sobolev approximation (see the Appendix). The following 

simplifications were used: 

1 . The scattering in the medium was assumed to be coher- 
ent (Thomson scattering). This approximationmay be 
considered justified, because the radiation incident on 
the disk has a fairly soft spectrum. The effect fromthe 
change in the frequency of photons due to the direct 
Compton effect as they are scattered in the disk will be 
small for photons with energy E < 20 keV (for more 
detail, see below). 

2. In solving the transfer equation, we use the Eddington 
approximation (see Eq. (A14)). 

3. We divided a plane-parallel layer of finite thickness 
into three zones in each of which the absorption co- 
efficient depends only on frequency: k = k{v). The 
central zone corresponds to a cold disk (n = k^), while 
the zones at the top and the bottom describe the atmo- 
spheric layer on both sides of the disk (k = Kam) (see 
Fig. 1). 

We assume that the opacity in the disk is determined by 
photoabsorption for a cold gas (Morrison and McCam- 
mon 1983): 

Kd = Kmm{J^)- (10) 



For the atmosphere, we will consider two cases: 

• a completely scattering atmosphere: 

nam{y) = (i- (11) 

• the absorption and scattering coefficients in the 
atmosphere are approximately equal: 



l^atm ~ 



f , Kmm > cr 



(12) 



where a is the scattering coefficient. 



Solving the transfer equation for X-ray photons by the 
Sobolev method, we obtained fairly simple analytic formu- 
las for the mean intensity and flux of X-ray photons at a 
given frequency as a function of the surface density (which 
is accumulated from the surface deep into the disk) S and 
parameters So, ^am, as well as the specified X-ray opaci- 



ties and K^„: 



(13) 
(14) 



where Sq is the total surface density in the disk. Here and 
below, the superscript v emphasizes the dependence of a 
quantity on frequency. 

The mean intensity and flux of X-ray photons in the disk 
at frequency u are (for more detail, see Appendix 1) 



+(1-0") 
//r„.(S)=F^-| — 



e-r7Co+g-(ro''-r'')/Co 



-r"ICo_g-(r^-r'')/Co 



(15) 



(16) 



where Tq is the total optical depth of the disk in the vertical 
direction for X-ray radiation at frequency v, = (a+ k^S, 
A = k'^ is the absorption coefficient for X-ray photons 
with frequency v in the disk, a is the scattering coefficient, 
/t= V'3(l-A). 

For a disk without an extended atmosphere, the coeffi- 
cients D'^ and C are defined by Eqs. (A23) and (A27), re- 
spectively. For a disk with an atmosphere, the coefficient 
D'^, as above, is calculated from Eq. (A23), while the co- 
efficient has a slightly more complex form (see (A49)). 
Note that D" is not an independent coefficient; it is related 
to C (see, e.g., Eq. (A27)). 

The additional heating of the disk by X-ray photons of a 
given frequency is proportional to their mean intensity: 



e'^=47rpK^/; 



1/ jl/ 

tot' 



(17) 
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Note that 7/^, tot includes both primary and scattered Xray 
photons. 

The local energy release in the disk through its irradiation 
by X-ray photons with a spectrum 8(1^) is 

POO POO 

e= e''du = 4np • (18) 

Ja Jo 

The total heating of the disk from the photosphere to the 
central plane through its irradiation is 

POO 

Qirr(^ph)= / H^JJ:pk)dl^, (19) 

Jo 

where E^/, is the surface density in the layers above the disk 
photosphere. 

The accuracy of the approximate solution of the transfer equation 
by the Sobolev method when calculating the held of X-ray 
photons in the disk. 

It is necessary to estimate the accuracy of the simple for- 
mulas for the flux and mean intensity of X-ray photons in the 
disk (17) and (18) derived by the Sobolev method. For this 
purpose, we compared the plane albedos of a semi-infinite 
layer Asob (Eq. (A39)) obtained by the Sobolev method for 
various incident X-ray photon energies (we used the energy 
dependence of the absorption coefficient Kmm(E) for a cold 
gas) with the exact albedos: 

• Aexact — the exact albedos for coherent scattering with 
the Rayleigh phase function were calculated via the 
Chandrasekhar //-functions (for more detail, see Ap- 
pendix 3). 

• '^sim'' — the albedos for scattering with the Comp- 
ton effect. To calculate A^,"™'', we used the method 
of Monte Carlo numerical simulations (see, e.g., Pozd- 
nyakov et al. 1983). The semi-infinite layer was as- 
sumed to be cold, kT <C w^c^; 10^ trial photons were 
taken for each albedo. To estimate the accuracy of 
our Monte Carlo simulations, we also determined the 
albedo by this method for coherent scattering A^m, 
which can be compared with its exact value of A^j^a^. 

Figure 2 shows the dependence of various albedos (Asob, 
^.lim ' °" energy E (upper plots) and ratios 
A sob/ -A exact, Asim/Aexact, and A^"Z'' j A^xact (lowcr plots) for 
two angles of incidence of the X-ray photons on the disk, 
cos^o = 0.1 (left) and cos^o =1-0 (right). We can estimate 
the accuracy of the derived approximate formulas (17) and 
(18) for the mean intensity and flux of X-ray photons in the 
disk as < 10% for photons with energy £ < 20 keV. 

Note that calculating the field of X-ray photons with a 
higher accuracy will require solving the problem of X-ray 
photon transfer in the disk by taking into account the fre- 
quency change in each scattering. This problem can be 
solved numerically (for a cold, geometrically thin disk, see, 
e.g., Psaltis 2002). 



23 The System of Equations for the Vertical Structure of 
a Disk Irradiated by a Central X-ray Source and its 
Solution 

The vertical structure of the accretion disk at a given ra- 
dius can be calculated under the assumption of a standard 
a-disk (Shakura and Sunyaev 1973); the (r,(^) component 
of the viscous stress tensor in the disk is proportional to the 
pressure = aP. Here, we assumed that a = 0.5. 

Denote the density, temperature, energy flux in the verti- 
cal direction toward the surface, surface density, (Rosseland) 
optical depth and opacity, and the photospheric level for in- 
trinsic radiation (from the disk midplane) by p, T, Q, S, r'', 
k'', and Zph, respectively. Let us introduce the vertical coor- 
dinate z' measured from the photospheric level deep into the 
disk (see Fig. 1). The surface density S is also measured 
toward the central plane of the accretion disk: 

E = Spft+/ pdz, (20) 
Jo 

where T,ph is the surface density in a column of material 
above the disk photosphere; the total surface density is Eq 
in the disk and E = Eo/2 in the its central plane. 

To solve the one-dimensional problem of the vertical ac- 
cretion disk structure at a given radius (R), we sought a nu- 
merical solution of the following system of ordinary differ- 
ential equations (see also Suleimanov et al. 2007): 

1. The equation of hydrostatic equiUbrium in z coordi- 
nate: 

^ = pfllizpb-z'). (21) 

The contribution from the radiation pressure is insignif- 
icant in the outer disk: Pgas ^ Prod- We neglect the 
contribution from Prad to the total pressure by assum- 
ing that P=Pjas. 

2. The equation of energy transfer by radiation in z coor- 
dinate (in the diffusion approximation): 

dT 3 k'' p 

Note that the Rosseland approximation (24) is appU- 
cable only for fairly deep layers, t'' > 1, under the 
disk photosphere. As long as the accretion disk is opti- 
cally thick in the vertical direction (tq ^1), using the 
diffusion approximation to determine the vertical disk 
structure is justified. At the same time, to find, for in- 
stance, the radiation spectrum from the accretion disk, 
the radiative transfer in the near-surface layers where 
the emergent radiation is formed must be calculated ac- 
curately (see, e.g., Suleimanov et al. 1999). 

It should also be noted that Eq. (24) disregards the 
energy transfer through convection. The convection 
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Fig. 2. Albedo of a semi-infinite layer with the opacity of a cold gas kmm versus energy E for two angles of incidence of the X-ray photons on the layer: 
cos 9 = 0.1 (left panels) and cos 6 = 1.0 (right panels). The upper plots show: Asot (solid line)-the approximate Sobolev value, Aexact (dashed line) — 
the exact value for coherent scattering, Aj,,,, (open circles) — from our numerical simulation for coherent scattering, A™""'' (filled circles) — from our 
numerical simulation with the Compton effect. The lower plots show the ratios A„i,/Ae-raci , Asj,,, /Ae_-raci snd A'^"!^'' /A„aa ■ 



mechanism is switched on when the vertical temper- 

dT 

ature gradient (24) is larger than the adiabatic — > 

dz 

dT 

— . The convection zone emerges in the outer ac- 

d?' 

ad 

cretion disk where the temperature drops below T ^ 
10* K and hydrogen transforms into a neutral state (see, 
e.g., Meyer and Meyer-Hofmeister 1982). A thermal 
instability develops in the zone with partial hydrogen 
ionization in the accretion disk. Since we consider only 
the case of stationary accretion, we calculated the disk 
structure up to the radius Rk at which the zone with 
incomplete hydrogen ionization appears. 



3. The equation for energy release: 



4. The equation for the surface density E: 



dQ 3 



(23) 



The sum of two terms appears on the right-hand side 
of Eq. (25). The first term represents the intrinsic en- 
ergy release in a standard Shakura-Sunyaev disk, while 
the second term represents the additional energy re- 
lease due to the absorption and thermalization of direct 
and scattered X-ray photons in the accretion disk (with- 
out irradiation e = 0). To find e dependent on depth z', 
we use an approximate solution of the X-ray radiation 
transfer problem (see Eq. (20)). 



dz' ^' 

5. We find the optical depth t'' from the equation: 

dr' 



(24) 



dz 



(25) 



The Rosseland opacity in Eqs. (24) and (27) depends 
on density and temperature: k'' = K''ip, T). For a given 
chemical composition, we determined n'' using the 
Opacity Project tables (Badnell et al. 2005). 

The equation of state (for an ideal gas) should be added 
to Eqs. (23)-(27): 

P = -pT. (26) 
p 

We determined the molecular weight p at the photospheric 
level of the disk using the Saha formula (for a given chemi- 
cal composition): 



p = piPphJph) 



(27) 



We found the radius Rk at which the zone with incomplete 
hydrogen ionization appeared in the accretion disk from the 
condition 



(28) 
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where /x,h = 0.66 and /z„; = 1 .26 are the molecular weights of 
the material in which only hydrogen is ionized (at the chem- 
ical composition specified above) and a completely neutral 
medium, respectively. We calculated the vertical structure 
of a stationary disk up to the radius Rk- 

To solve the system of ordinary differential equations 
(23)-(29) with (28) and (29), we specified the following 
boundary conditions at the photospheric level of the disk: 



where pph = 



Tph 

Qph 

^ph 
fJ-Pph^ 

^Tph 



X T^I, + Pnl 



^K^ph 
K^PphJph)^''"' 

QphV" 



o'sb) 

+ Qirri^ph) 



(29) 



'ph 



K'iPphJph) 

2 



Paim, Qvis, Qirr arc, respectively. 



the surface density and pressure of the atmospheric layer 
(see Eqs. (8) and (10)), the total energy release in the disk 
through viscosity, and the total energy release through X-ray 
heating of the disk; asB is the Stefan-Boltzmann constant. 
We derived the dependence 2,rr(S) from the solution of the 
X-ray radiation transfer problem (see (21)). We have 2,>r = 
0, J^atm = 0, for an unirradiated disk and Qirr > 0, T,atm = 
for an irradiated disk without an atmosphere. 

The intrinsic energy release in the accretion disk at radius 
R is defined by the expression 



Qv 



3GMM 
87r/?3 



-f(R) 



(30) 



where the function f(R) is related to the boundary condi- 
tion at the inner accretion disk boundary. We investigate the 
outer disk where it can be assumed that /(R) = 1. 

Given the boundary conditions (31), we integrated the 
system of equations (23)-(27) from the photospheric level 
iz' = 0) deep into the disk surface using the Runge-Kutta 
method implemented in the Numerical Recipes software 
package (Press et al. 1992). The height of the photosphere 
Zph is a free parameter of the problem and is not known in 
advance. For an irradiated disk, the total surface density in 
the disk Eo is another free parameter. We found Zph and So 
when constructing the solution by the method of successive 
iterations using an additional condition for the energy flux 
in the central plane of the disk: 



Q 



■ 0. 



(31) 



z'=Zph 



3 RESULTS AND DISCUSSION 

Figures 3-6 present the vertical structure of an irradi- 
ated accretion disk (the dependences Qiz), T(z), P{z), E(z), 



T^z) and k''{z)) for two radii, R = 10^" cm and 10^1 cm, and 
two luminosities of the central Xray source, = 0.1 and 
,1.0LEdd- The plots show the foUowing models. 

1. the irradiation of a disk with a completely scattering 
atmosphere {Kam = 0) (thick solid line); 

2. the irradiation of a disk with an atmosphere where ab- 
sorption is approximately equal to scattering for pho- 
tons with energy E < lOkeV (see Eq. (14); thin soUd 
Une); 

3. the irradiation of a disk without an atmosphere (dashed 
fine); 

4. an unirradiated disk (dotted Une). 

We considered only the models of a stationary accretion 
disk, i.e., those satisfying the following condition: the tem- 
perature on the disk photosphere should be higher than the 
hydrogen recombination limit. For this reason, we did not 
provide the model of an unirradiated disk for the radius 
R = 10'^ cm and the calculation of an irradiated disk with 
an atmosphere is shown only for Lx = ifidd- On the plots 
of Qiz) atR= 10'^ cm, instead of model (4) the dotted line 
indicates the vertical profile of the energy flux due to vis- 
cous heating in the disk, Qyisiz), for a completely scattering 
atmosphere. 

Note that the heating of the accretion disk through irra- 
diation at R = 10^" cm affects up to half of its height. At 
R= 10'^ cm, irradiation can have an effect on the entire ver- 
tical disk structure, up to the central plane. The accretion 
disk is deeply heated by the scattered X-ray photons, while 
the direct photons from the central source are absorbed and 
thermaUzed in a thin near-surface layer. The heating effect 
is enhanced significantly in the presence of an extended at- 
mospheric layer above the accretion disk. 

Mostly hard X-ray photons, which can traverse a suffi- 
cient distance in z coordinate without absorption, are in- 
volved in heating the deep layers of the accretion disk. In- 
deed, as follows from Eq. (17), the mean intensity of X-ray 
photons at a sufficiently large depth S from the surface of 
a semi-infinite layer is determined by the field of scattered 
photons: 



7(£;S) cx C(£)exp(-S/Srf,/(£)) 



(32) 



Figure 7b shows the dependence of the quantity C on X-ray 
photon energy E and, for comparison, the spectrum of the 
radiation incident on the disk S{E). Figure 7a shows the 
characteristic surface density E^^,/ (the intensity of scattered 
radiation with photon energy E at this depth decreases by a 
factor of e; see (34)) as a function of E. The quantity S^,/ 
was calculated for the characteristic angle of incidence of 
the X-ray photons cos^o = 0.1 using an energy dependence 
of the absorption coefficient like that for a cold gas (Morri- 
son and McCammon 1983). 
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Fig. 3. Vertical structure of an irradiated accretion disk: R = lO'" cm, Lx = Z-Edd- The central object is a neutron star with mass Mi = IAMq, the viscosity 
parameter in the disk is a = 0.5. The thick solid line indicates the irradiation of a disk with a completely scattering atmosphere, the thin solid line indicates 
the irradiation of a disk with an atmosphere where absorption is equal to scattering for photons with energy £ < 10 keV, the dashed Une indicates the 
irradiation of a disk without an atmosphere, and the dotted line indicates an unirradiated disk. 
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Fig. 4. Vertical structure of an irradiated disk: R = lO"* cm, Lx = O.liEdd- Tlie designations of tlie lines for variousmodels of an irradiated disk are the 
same as those in Fig. 3. 



ASTRONOMY LETTERS Vol.37 No. 5 2011 



10 



Mescheryakov et al. 




Fig. S. Vertical structure of an irradiated dislc: R= 1 " 6l, Lx = Ltdd ■ Tiie designations of the lines for variousmodels of an irradiated disk are the same 
as those in Fig. 3. The model of an unirradiated disk is not shown. Instead, the dotted line on the plot of Q(z) indicates the profile of the energy flux due to 
viscous heating (Qvij) in the model of an irradiated disk (a completely scattering atmosphere). 
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Fig. 6. Vertical structure of an irradiated dislc: R= lO" 6l, Lx = 0. 1 LEdd- The designations of the Unes for various models of an irradiated disk are the 
same as those in Fig. 3. The models of an unirradiated disk and an irradiated disk without an atmosphere are not shown. The dotted line on the plot of Q(z) 
indicates the profile of the energy flux due to viscous disk heating in the model of a disk with an atmosphere. 
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radius where hydrogen recombination begins (and, ac- 
cordingly, Zph decreases sharply). 

• In the case of disk irradiation, the zone with incomplete 
hydrogen ionization is greatly shifted toward large 
radii: Rk increases by a factor of ~ 10 in an irradiated 
disk (with an atmosphere) compared to an unirradiated 
one. 




1 2 3456789 10 20 

6 J — I 1 1 — I — I — r~rn 1 




-4 I I I I I I I I I 

1 2 3456789 10 20 

E , keV 



Fig. 7. Energy dependence of the characteristic penetration depth of scat- 
tered X-ray photons Srf,y(£) (top panel) and the quantity C(E) (bottom 
panel) (a semi-infinite layer, cos$ = 0.1). For comparison, the dotted Une 
indicates the incident X-ray spectrum. 

It follows from Fig. 7 that the mean intensity of scat- 
tered X-ray photons with energy E=10 20 keV is attenu- 
ated by a factor of e at a depth S^,/ « 1 4 g cm"2 (which 
corresponds to an optical depth in the disk for intrinsic ra- 
diation r'' « 10 -=-40). Therefore, the hard X-ray radiation 
from the central source (with energy E >1Q keV for a cold 
disk) at the outer disk radius, where the total surface density 
in the vertical direction reaches Sq ~ 20 g cm"2, is capable 
of heating all disk layers up to the central plane (see Figs. 
5 and 6). At the same time, the accretion disk at this radius 
still remains optically thick, t'' > 100, for intrinsic radiation. 

In Fig. 8, the quantities Q/Qvis, Tc, Tph and Zph (the ra- 
tio of the total flux to the viscous one at the photospheric 
level, the temperature in the central plane of the disk and on 
the photosphere, and the height of the photosphere, respec- 
tively) are plotted against radius R for various models of an 
irradiated disk. The radial profiles are shown up to the ra- 
dius at which a zone with incomplete hydrogen ionization 
appears in the disk. 

We will note the following: 

• the difference between the temperatures in the central 
plane and on the photosphere of an irradiated accre- 
tion disk decreases appreciably with increasing radius: 
from Tc/Tph « 4.5 at ^ = 10'* cm to Tc/Tph « 1.5 at 
/? = 10" of. At /? > 10" cm, the vertical structure of 
an optically thick disk is essentially isothermal. 

• Even strong irradiation has no significant effect on the 
disk thickness at the photospheric level Zph up to the 



• For the luminosity of the central source Lx = Lndd. the 
outer radius of a stationary irradiated disk (with an at- 
mosphere) can be « 6 X 10" cm, in qualitative agree- 
ment with the estimates of the outer radius in long- 
period persistent LMXBs. 

• In the model of an irradiated disk with an atmosphere, 
the problem of heating the outer (R > Rg) cold accre- 
tion disk that is screened (see Dubus et al. 1999) from 
the direct photons of the central Xray source is natu- 
rally resolved. 

4 CONCLUSIONS 

Here, we investigated the vertical structure of the outer ac- 
cretion disk in LMXBs by taking into account the possibility 
of the scattering of X-ray photons from the central source in 
the disk and the atmospheric layer The atmospheric param- 
eters were taken from Jimenez-Garate et al. (2002). We 
reached the following conclusions. 

1. We derived simple analytic expressions for the disk 
heating by scattered X-ray photons using an approxi- 
mate solution of the transfer equation by the Sobolev 
method. This approximation has a < 10% accuracy in 
the range of X-ray photon energies E < 20keV. 

2. We showed that the scattering of X-ray photons by free 
electrons affects significantly the heating of the outer 
accretion disk in LMXBs. Having been scattered, the 
X-ray photons with energy E > lOkeV incident from 
the central source at a small angle to the disk surface 
can penetrate fairly deep into the disk photosphere and 
can affect significantly the vertical structure of the ac- 
cretion disk at outer radii. 

3. The scattering of X-ray photons is particularly impor- 
tant at large radii, where the total disk surface density is 
fairly low Sq < 20 g cm~^ (in this case, the disk still re- 
mains optically thick, r'' > 100). At radii R> lO" cm, 
irradiation can heat all layers of an optically thick disk 
and its vertical structure becomes essentially isother- 
mal. 

4. When a disk with an atmosphere is irradiated, the ra- 
dius at which a zone with incomplete hydrogen ion- 
ization appears in the disk increases by a factor of 
~ 10 compared to an unirradiated disk. For exam- 
ple, for the luminosity of the central source Lx = L^^, 
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Fig. 8. Radial structure of an irradiated accretion disk for two luminosities of tlie central X-ray source: Lx = Z-Edd (left panels) and Lx = O.lLEdd (right 
panels). The upper row of plots shows the ratio of the total energy flux in the vertical direction to the energy flux only due to viscous heating Q/Qris at the 
photospheric level as a function of radius R: the middle rowsh ows the temperature in the central plane Tc and on the photosphere Tpi, of the disk; the lower 
row shows the height of the photosphere Zph/R and atmosphere Zaim/R (taken from Jimenez-Garate et al, 2002; the hne with long dashed on the plots). 
The designations of the lines for various models of an irradiated disk are the same as those in Fig. 3. 



the outer radius of a stationary irradiated disk can be 
/? ss 6 X 10^^ cm, in agreement with the estimates of the 
outer disk radius in the long-period persistent LMXBs 
GX 13H-1 and CygX-2. 

5. In the model of an irradiated accretion disk with an at- 
mosphere, the problem of the heating of its outer cold 
(T < 10"^ K) parts, which are screened (see Dubus et 
al. 1999) from the direct photons of the central X-ray 
source, is naturally resolved. 
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A ANALYTIC SOLUTION OF THE PROBLEM OF 
X-RAY PHOTON TRANSFER IN A 
PLANE-PARALLEL LAYER IN THE SOBOLEV 
APPROXIMATION 

A.l A Layer with a Constant Absorption Coefficient ( a 
Disk without an Atmosphere) 

Consider a plane-parallel layer of material with a constant 
(in depth) absorption coefficient k on which a parallel X- 
ray flux is incident at an angle to the inward normal to 
the surface. The same flux is symmetrically incident on the 
lower surface at an angle (180° -^o)- The azimuth of the 
incident radiation is zero, (f)o = 0. 

Consider the irradiation of a cold layer with temperature 
T by X-ray photons with frequency ly, so that the inequality 
kT <^ hu is vaUd; there are no intrinsic energy sources in 
the layer We assume the scattering in the medium to be 
coherent (Thomson scattering). In this case, the intensity 
of scattered radiation at frequency v, I{t, C, (!>), satisfies the 
integro-differential equation (see, e.g., Chandrasekhar 1950) 

C^^^^^=-Iir,(:,4>)+SiT,C,cf>), (91) 
The source function S(t, (fi) is 

S(T,C,(p) = Sinc(T,C,<P)+ 

(52) 

where r is the optical depth in the vertical direction mea- 
sured from the surface deep into the layer, ( = cos 9, X= 
is the single-scattering albedo, k is the absorption coeffi- 
cient at frequency v, a is the scattering coefficient, and 
x((,<p,(' t'P') is the phase function. 

In the case of randomly oriented scatterings, the phase 
function is a function of the scattering angle alone, x(cos 7), 

cos7 = CC' + Vl-CVl-C'^cos((/)-^')- (93) 

We consider the scattering by free electrons with the phase 
function 



x(cos7) =4(1 + cos 7 ) = 1 +X2/'2(cos7), (94) 

where F2(cos7) = ^"2''' the second Legendre polyno- 
mial. The coefficient X2 = 1/2 for the Rayleigh phase func- 
tion and X2 = for isotropic scattering. 

The term 5,„c in the source function (A2) allows for the 
single-scattering of the direct radiation incident on the layer: 

Sinc(T, C> 0) = ^ / / -^Ccos 7)/;„c(t, C , (i)')dC'd(j)' , 

(95) 

where /,„c(r, (,(p) is the intensity of the direct radiation in 
the layer:: 



Co = cos 6*0, S(x) is the Dirac delta function, and Fx is the 
monochromatic iUumination of the area perpendicular to the 
direction of incidence of the X rays. Thus, 



Sinc(T,C,<l>) = ^•«(COS70) 



g-'r/Co_|_g-(To-r)/Co 



where 



cos 7o = CCo + Vl-C^V'l-Co cos . 



(97) 



(98) 



For an approximate solution of the transfer equation for 
scattered photons (Al), (A2), (A7), we use the method de- 
veloped by V. V. Sobolev for the problem of diffuse reflec- 
tion and transmission (Sobolev 1956). This approximate 
method allows fairly simple formulas to be derived for the 
mean intensity of X-ray photons in a layer J and the albedo 
A (see below and the textbook by Nagirner 2001). 

The intensity I{T,Q,(j)) and the source function S{T,(,(t>) 
can be written as a decomposition into azimuthal harmonics. 
For the Rayleigh phase function, we have 

I{T,C,<t>)= k{T, + 2/1 (T, COS <j> + 2/2(r , COS 2<p, (99) 



S{t, C, (p) = SoiT, + 2Si{t, cos 0+ 252(r , Q cos 2<l>. 

(910) 

The mean intensity and flux of scattered X-ray photons in 
the vertical direction are defined by the zeroth azimuthal har- 
monic of the intensity Iq(t, Q: 



1 



Ioir,OdC, 



(911) 



/+1 nlTT 
^ C,I{TX,4»d(:d<t> 

= 27r£ CIoiT,OdC. (912) 

The transfer equation for /o(t, Q is 

C^^^P^=-Ur,0+So(T,0. (913) 

OT 

We use the Eddington approximation 

1 



2 7-1 



1 



MEdd 



(914) 



+Fxe 



-(To-T)/Co 



'5(C + Co)<5(<^), (96) 



where /XEdd is the constant Eddington factor. In the general 
(more accurate) case, /XEdd changes with z coordinate. 
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Substituting (A3), (A4), (A7), and (A9) into (A2), us- 
ing the Eddington approximation (A14), and gathering the 
azimuth-independent terms, we have for the source function 



5o(t,C) = XJ+ 



l + ^(l-3Co^)(l-3a 



4n 

T"/Co + g-(-ro-r)/Co 



(915) 



X2 = 1/2 in the case of scattering with the Rayleigh phase 
function and X2 = in the case of isotropic scattering. 

Let us integrate the transfer equation for the zeroth har- 
monic (A 13) over ( with the weight factors 1 and i^. Apply- 
ing again the Eddington approximation, we will obtain the 
following system of differential equations for the functions 
Jandh=^: 



The coefficient/) can be found by directly substituting solu- 
tion (A22) into Eq. (A18): 



D = 



(523) 



The general solution of the homogeneous equation (A21) 

is 

y,(.)=^£x(^'^"";^^^" ■^</''^>0 (524) 
^ 4n \ C1T + C2 ,A= 1,«; = 

We are interested in the case of k > 0, A < 1. 

Since the problem is symmetric relative to the r = ro/2 
plane, we will obtain 



JJt): 



CFx 
4tt 



(525) 



^=_(l_A)y.^ 
dr 4tt 



/Co + g-(TO-T)/Co 



(516) 



To find the coefficient C, we use the condition at the outer 
boundary of the layer r = 0: 



h--- — 
3d^ ■ 



(617) 



Substituting h from (A17) into (A16), we will obtain an 
inhomogeneous second-order differential equation for the 
monochromatic mean intensity of scattered photons: 



■k^J = 



3XFx 
4tt 



g-T/C0_|_g-(T-0-T)/C0 



(518) 



where k = \/3(l-A). Having determined J(t), we can find 
the flux of scattered photons //(r) from Eq. (A17): 



H-- 



4-KdJ 
'Ydr' 



(519) 



Interestingly, Eqs. (A18) and (A19) derived by the 
Sobolev method have the same form for the isotropic or 
Rayleigh phase function (the parameter X2 does not enter 
into the formulas). 

A general solution of the inhomogeneous differential 
equation (A18) should be sought in the form of the sum 



JiT) = JpiT)+JgiT), 



(520) 



where Jpir) is a particular solution of the inhomogeneous 
equation (A 18) and Jgir) is the general solution of the ho- 
mogeneous equation 



dr^ 



(521) 



The particular solution of the inhomogeneous differential 
equation (A18) is 



Jp = - 



4tt 



'"/Co _|_ g-(To-r)/Co 



(522) 



J- 



2dJ 
3d^ 



= 0. 



(526) 



T=0 



This is equivalent to H\^^q = -2ttJ\^^q — there is no ex- 
ternally incident scattered radiation at the boundary of the 
layer. Substituting (A20) into condition (A26) yields an ex- 
pression for the coefficient C: 



C = Dx 



1 + g-^o/ Co + ( 1 + g-T-o/ Co 

3 Co 



(527) 



Thus, the mean intensity and flux of scattered Xray pho- 
tons in the layer at an optical depth r are 



J = 



Fx 
47r 

Fx 
4n 



g-T"/Co+g-(T-0-T)/Co 



(528) 



H = 



kC 

T 

D 

3Co 



Fx 

Fx^ 



g-kr _ g-k(To-T) 



g-T/Co _ g-(-ro-r)/Co 



(529) 



where the coefficients D and C are determined from Eqs. 
(A23) and (A27), respectively. 

We emphasize that D is not an independent coefficient; it 
is related to the coefficient C (see (A27)). We are dealing 
with a second-order differential equation that has two inde- 
pendent solutions. For the semiinfinite problem, the second 
solution is ruled out, because it grows exponentiaUy. For a 
layer with a finite optical depth (tq), we naturally have a so- 
lution with two constants Ci and C2 that, in our case, turn 
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into one constant C, because the upper and lower parts of 
the layer are irradiated symmetrically. 

The above expressions for the mean intensity J(t) and 
flux H(t) of scattered photons in the layer are so structured 
that we have an indeterminate form of the type (oo - oo) at 

A. = l-T^, (530) 

because the denominator of the coefficients D and C be- 
comes zero. Let us evaluate this indeterminate form. At 
A = A*, the functions 7(t) and H{t) will then have finite val- 
ues: 



y(A.;T) = -^(C-Z)),.i 

47r Co 



/?(A.;r)=^(C-£»)^^ 



e r/Co+g-(r„-r)/Co I (53 



""/Co _ ^-(to-t)/Co 



(532) 



3 Co 



1 + 



3to 
2 



1 + -^ 

3Co 



-To/Co 



1+g-To/Co 



(533) 



If the X-ray photons are incident at a fairly small angle 
(90° -^o) to the layer surface, more specifically, (90° -6*0) < 
arcos{\/ \/3) w 35.3°, then the condition kC,Q < 1 is always 
met and the denominator of the coefficient D does not be- 
come zero for any A. 

It is important to note that there must also be an angle 
of incidence in the exact solution of the transfer problem 
with external irradiation where a similar indeterminate form 
emerges. It stems from the fact that at large angles of inci- 
dence the field of scattered photons decreases with depth 
more rapidly than the field of incident photons, while at 
small angles the field of incident photons decreases more 
rapidly. 

Thus, the Eddington approximation (/iEdd = 1/3, see 
(A14)) gives a quahtatively newresu It (the appearance of 
an indeterminate form at some critical angle of incidence) 
in problems with external irradiation. As the critical angle 
is approached, the constants C and D tend to infinity, but 
their difference C-D remains finite. Therein lies the deep 
meaning of the approximate (Sobolev) solution. This result 
must also be retained in the exact solution of the problem. 

The mean intensity and flux in the vertical direction of the 
direct radiation incident on the layer are 



(534) 
(535) 



The total mean intensity and flux of X-ray photons in the 
layer at an optical depth r (direct and scattered radiation) 
are 



Jtnt — ~. C 



4lT 

Fx 
4tt 



g-kr + g-kiro-T) 



+ 



(l-D) 



(536) 



= F. 



kC 



+ 

e ^ —e Co 



(537) 



If the X-ray photons are incident at a small angle to the sur- 
face, then the direct radiation does not pass deep into the 
layer, because the exponential e"^/'-", 1/Co ^ 1, decreases 
rapidly and the radiation field inside the layer is determined 
by diffuse radiation. 

The albedo of the layer for X-ray photons is 



H 



(638) 



T=0 



Substituting H and from (A29) and (A35) yields 



A = 



D kC f \ 



-kro 



3Co Ko\l-e-^»/'^o 



(639) 



To estimate the accuracy of the derived formulas for the- 
mean intensity and flux of X-ray photons (A36) and (A37), 
we compared the albedos (A39) with their exact values for 
various angles of incidence of the Xray photons in the case 
of a semi-infinite medium (tq 00). The exact albedos for 
the isotropic and Rayleigh phase functions were calculated 
via the Chandrasekhar //-functions (for more detail, see Ap- 
pendix 3). 

The table gives approximate albedos obtained by the 
Sobolev method for various ( and A, and their exact val- 
ues for isotropic and Rayleigh scattering. We see excellent 
agreement between the exact and approximate albedos. 

A.2 A Plane -Parallel Layer with Two Absorption 
Coefficients (Disk+ Atmosphere) 

Let us divide a plane layer into three zones at the bound- 
ary of which the absorption coefficient changes abruptly (see 
Fig. 1): K = Ka 'ivi zones 1 and 3 and k = in zone 2. The 
middle zone 2 corresponds to a cold disk, while zones 1 and 
3 correspond to atmospheric regions at the top and the bot- 
tom. The central plane of the disk passes in the middle of 
zone 2 and is the plane of symmetry of the problem. Denote 
the optical depth in the vertical direction in zones 1 and 3 by 
Tfl and the optical depth between the boundary of zone 2 and 
the plane of symmetry of the layer by Td (see Fig. 1); the 
total optical depth in the layer is tq = 2{Ta + Td). 

Inside each zone, we seek a solution of Eq. (A18) with a 
constant A: 



\a ,0<T<Ta 
X=( Ad ,Ta<T <To-Ta 
K ,TQ-Ta<T<TQ. 



(540) 



where Xa = Xd = 

It will suffice to find a solution in the region < r < to/2. 
Since the problem is syrmnetric relative to the t = tq/2 
plane, the solution for the symmetric zone to/2 < r < tq 
can be obtained by the change of variables r <;=^ tq - r. 
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For XaondXd < 1 iKa,Kii> 0), the general solution of the 
inhomogeneous differential equation (A18) in zone 1, < 

r < Ta. 



J{t) 



El. 



——Da \ e '^0 +e «o 



(541) 



ka = V3(1-AJ. 

In zone 2, r,, < r < to/2, given the symmetry relative to 
the ro/2 plane, we have 



J{t) 



F,Cd ( 



An 



— - — [ e ''0 +e ''0 
An \ 



(542) 



fcrf = V3(l-Ad). 

The coefficients Da and D^ are defined by the expressions 



Da = 

Dd = 



3 Co 
^-k'dCl 



(643) 
(544) 



To find the coefficients Cia, C2a, and Q, it is necessary to 
use the following three boundary conditions: 

• the condition at the upper boundary of the layer (A26); 

• at the boundary of zones 1 and 2, the quantities k, A and 
k undergo a discontinuity, while the quantities J and H 

must change continuously. Hence two joining condi- 
tions for the mean intensity and its derivative follow: 



J 

dJ 

J? 



= J 



T=Ta-Q 



T=Ta-0 



dJ 



dT 



(945) 
(546) 



T=r„+0 



Substituting (A41) and (A42) into conditions (A26), 
(A45), and (A46), we obtain a system of three linear alge- 
braic equations with three unknowns solving which we find 

Cia, Cla, and Q: 



Daf 



ka 



e-kar.-iPa-Dd) 



g-b+ _ b- 
g+ Co 



1 -^k 



l + jka 



a. ""fl 



'-+ka 

g+ " 



Cia ■ 



Daf ^ 1-ffea 



l + lka 



l + lka 



Cd- 



Cue-'""^' +C2ae*''^° - {Da -Dd)K 



(547) 

(648) 
(649) 



The auxihary quantities /, b+, b-, g+ and g- are defined by 
the formulas: 



/ = 

K = 
b- = 



2 



— Z2. 

— la. 



g-kdTa _|_ ^-kd{To-Ta) 



e ^ 



Co 



Co 



g- = kdx (e-'''^'-e-'''^^''-^''>y 

As above, the flux of scattered photons H(t) and the 
albedo of the layer A can be found using Eqs. (A 19) and 
(A38). 

A.3 The Albedo of a Semi-infinite Layer: Exact Values 

The albedo of a semi-infinite atmosphere in the case of co- 
herent scattering can be calculated using the Chandrasekhar 
//-functions (Chandrasekhar 1950). The //-function can be 
found as a solution of the nonlinear integral equation 

i){ix)H{ri)d,i 



h(m)=i+m^^(m) 



/i+77 



(650) 



where is the characteristic function (Chandrasekhar 
1950) for the chosen phase function x(/i); ^ ^ 1. Con- 
sider two cases: (1) isotropic scattering and (2) Rayleigh 
scattering. 

In these cases, the characteristic function is (Sobolev 
1968) 



A 



l + ^(3(l-A)/i' 



■1)/'2(M) 



(351) 



Piin) is the second Legendre polynomial, X2 = for isotropic 
scattering, X2 = 1/2 for Rayleigh scattering. The character- 
istic function (A51) satisfies the condition J^ip{r])dr] < i. 

We found the function //(/i) using the procedure of suc- 
cessive iterations proposed by Bosma and Rooij (1983) 
(see method 3). Having determined //(/x), we found the 
albedo of a semi-infinite atmosphere Ap(Q from the formu- 
las (Sobolev 1968) 



1-A 



ApiO -- 


= 1- A 


A = 


= M1N2-I 


Ml = 


- 1--/ 
2 Jo 


M2 = 


A 



H(0{N2-NiO, 



(952) 



Ni 
N2 



^X23il 



Hiv) 
Hin) 

Jo 



df] 



TjdT] , 

H{r])P2{v)vdv , 



= M, 



The albedos of a semi-infinite atmosphere for various Co and 
A for the isotropic and Rayleigh phase functions are given in 
the table. 
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Table. Albedos for a semi-infinite layer 
in the Sobolev approximation (As) and tiieir exact values for the isotropic (A,jo) and Rayleigh (A;,™) phase functions. 





C = 0.01 


C = 0.05 


C = o.io 


C = 0.20 


C = 0.40 


C = 0.60 


C = 0.80 


C = 1.00 


A 


As Ai,o Arl 




As Aiso Ah 


As Aiso Arl 






As Aiso Arl 





0.10 0.047 0.049 0.049 0.044 0.044 0.043 0.041 0.040 0.039 0.036 0.034 0.033 0.029 0.026 0.026 0.024 0.022 0.022 0.021 0.019 0.019 0.018 0.016 0.018 

0.20 0.097 0.101 0.101 0.091 0.091 0.090 0.085 0.083 0.081 0.075 0.071 0.069 0.061 0.056 0.055 0.051 0.047 0.047 0.044 0.040 0.041 0.039 0.035 0.038 

0.30 0.150 0.157 0.157 0.142 0.143 0.142 0.133 0.130 0.128 0.118 0.112 0.110 0.097 0.090 0.088 0.082 0.075 0.075 0.071 0.065 0.067 0.062 0.057 0.061 

0.40 0.208 0.218 0.218 0.198 0.199 0.198 0.186 0.183 0.181 0.166 0.159 0.157 0.137 0.128 0.127 0.117 0.109 0.109 0.102 0.094 0.096 0.090 0.083 0.088 

0.50 0.272 0.284 0.284 0.259 0.262 0.261 0.245 0.242 0.240 0.221 0.213 0.210 0.185 0.174 0.172 0.159 0.148 0.148 0.139 0.130 0.132 0.124 0.115 0.121 

0.60 0.343 0.358 0.358 0.329 0.333 0.332 0.313 0.310 0.308 0.284 0.276 0.273 0.241 0.229 0.227 0.209 0.197 0.197 0.185 0.174 0.176 0.165 0.155 0.161 

0.70 0.425 0.442 0.443 0.409 0.415 0.415 0.392 0.390 0.389 0.360 0.352 0.350 0.311 0.299 0.297 0.273 0.261 0.260 0.244 0.232 0.234 0.220 0.209 0.214 

0.80 0.524 0.543 0.545 0.508 0.516 0.517 0.490 0.4910.490 0.457 0.451 0.449 0.403 0.391 0.389 0.360 0.348 0.347 0.326 0.313 0.315 0.297 0.285 0.290 

0.90 0.656 0.675 0.678 0.642 0.652 0.654 0.625 0.629 0.630 0.594 0.592 0.591 0.541 0.532 0.531 0.496 0.486 0.485 0.458 0.447 0.449 0.426 0.415 0.418 

0.95 0.752 0.770 0.773 0.741 0.751 0.754 0.727 0.733 0.734 0.701 0.701 0.701 0.654 0.649 0.648 0.613 0.605 0.605 0.576 0.568 0.569 0.544 0.536 0.537 

0.99 0.886 0.897 0.899 0.880 0.887 0.889 0.872 0.878 0.879 0.858 0.860 0.861 0.830 0.829 0.830 0.804 0.802 0.802 0.780 0.776 0.776 0.756 0.753 0.753 



